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Abstract 

We analyze the octet baryon masses and the pion-nucleon (T-term in the framework 
of heavy baryon chiral perturbation theory. In contrast to previous investigations, 
we include all terms up-to-and-including quadratic order in the light quark masses. 
The pertinent low-energy constants are fixed from resonance exchange. This leaves 
as the only free parameter the baryon mass in the chiral limit, m. We find m= 
749ibl25 MeV together with (T7rAr(0) = 48±10 MeV. We discuss various implications 
of these results. 



1 Introduction 



The analysis of the octet baryon masses in the framework of chiral perturba- 
tion theory already has a long history, see e.g. [1-12]. In this paper, we present 
the results of a first calculation including a// terms of order (9(m^), where rrtq 
is a generic symbol for any one of the light quark masses mu^d,s- We work in the 
isospin limit m„ = irid and neglect the electromagnetic corrections. Previous 
investigations only considered mostly the so-called computable corrections of 
order or included some of the finite terms at this order. This, however, 
contradicts the spirit of chiral perturbation theory (CHPT) in that all terms 
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at a given order have to be retained, see e.g. [13-15]. In general, the quark 
mass expansion of the octet baryon masses takes the form 

m^m + Y.Bqm^ + Y.'^i^T + Y.Di^l + ■■■ (1) 
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modulo logs. Here, m is the mass in the chiral limit of vanishing quark masses 
and the coefficients Bq,Cq, Dq are state-dependent. Furthermore, they include 
contributions proprotional to some low-cncrgy constants which appear beyond 
leading order in the effective Lagrangian. In this letter, we evaluate these co- 
efficients for the octet baryons A, E and S. In addition, we also calculate 
the pion-nucleon cr-term which is intimately related to the quark mass expan- 
sion of the baryon masses [3,7,16]. For some comprehensive reviews, see e.g. 
[17-20]. 



2 Effective Lagrangian 



To perform the calculations, we make use of the effective meson-baryon La- 
grangian. Our notation is identical to the one used in [7] and we discuss here 
only the new terms. Denoting by the pseudoscalar Goldstone fields (tt, K, rf) 
and by B the baryon octet, the effective Lagrangian takes the form 



where the chiral dimension (i) counts the number of derivatives and/or meson 
mass insertions. The baryons are treated in the extreme non-relativistic limit 
[21,22], i.e. they are characterized by a four-velocity f^. In this approach, there 
is a one-to-one correspondence between the expansion in small momenta and 
quark masses and the expansion in Goldstone boson loops, i.e. a consistent 
power counting scheme emerges. The form of the lowest order meson-baryon 
Lagrangian is standard, see e.g. [7], and the meson Lagrangian is given in 
[15]. The dimension two meson-baryon Lagrangian can be written as (we only 
enumerate the terms which contribute) 
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C%-L^^r'^ + Y.hOf\ (3) 



with the 0\ monomials in the fields of chiral dimension two. Typical forms are 
Tr(i?[M^, [u^', 5]]), Ti{B[v ■ m, [w ■ m, B]]) or BBTT{Uf,u''), with = iu^d^Uu\ 
u — \/TJ and U — exp{i(f)/Fp) collects the pseudoscalars. The form of 
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is [7], 

^(2.tand) ^ Tr(5{x+, B}) + 6^ Tr(5[x+, B]) + ho Tr{BB) Tr(x+) , (4) 

i.e. it contains three low-energy constants and x+ = 'U^X'U^ + ux^u is propor- 
tional to the quark mass matrix M. — diag(mu, nid, mg) since % = 2BM.. Here, 
B — — (0|^g|0)/Fp and Fp is the pseudoscalar decay constant. All low-energy 
constants in d^^\ are finite. A subset of the hi has been estimated in [23] by 
analyzing kaon-nucleon scattering data. The splitting of the dimension two 
meson-baryon Lagrangian in Eq.(3) is motivated by the fact that while the 
first three terms appear in tree and loop graphs, the latter ten only come in 
via loops. There are seven terms contributing at dimension four, 

c%^j:diOt\ (5) 

i=l 

with typical forms of the O^^"* are BBTi{x\) Tr(i?[x_|_, At this 

stage, we take m„ = ^ m^. For m„ ^ ma, there is an additional term 
in jC^^^- The explicit expressions for the various terms in Eqs.(3,5) can be 
found in [24]. We point out that there are 20 a priori unkown constants. In 
addition, there arc the F and D coupling constants (subject to the constraint 
F + D = qa = 1.25) from the lowest order Lagrangian C^^^- What we have to 
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calculate are all one-loop graphs with insertions from C)^^ and tree graphs 

from C^^B^ ■ We stress that wc do not include the spin-3/2 decuplet in the 
effective field theory [6], but rather use these fields to estimate the pertinent 
low-energy constants (resonance saturation principle). We therefore strictly 
count in small quark masses and external momenta with no recourse to large 
Nc arguments. 



3 Baryon masses and pion— nucleon cr— term 

The form of the terms ~ rUg and ~ m^/^ for the baryon masses and 0-^^(0) 
is standard, we use here the same notation as Ref.[7]. The q'^ contribution to 
any octet baryon mass uib takes the form 

+ 4,B Ml, In(^) + Ml In(^) , (6) 

with P,Q = TT, K,ri and A the scale of dimensional regularization. The explicit 
form of the state-dependent prefactors ej,s can be found in Ref.[24]. Notice 
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the appearance of mixed terms ~ Mp Mq which were not considered in most 
existing investigations. The fourth order contribution to the baryon masses 
contains divergences proportional to (using dimensional regularization) 

with 'Je = 0.572215. These are renormalized by an appropriate choice of the 
low-energy constants di, 

di^d\{\) + V,L . (8) 



In fact, all seven di are divergent. The appearance of these divergences is in 
marked contrast to the calculation which is completely finite (in the heavy 
fermion approach). In what follows, we set A = 1 GeV. Similarly, the fourth 
order contribution to the pion-nucleon cr-term can be written as 



411(0) 



Ml 



A2 



3,<T 



ln( 



Ml ^ 
A2 



(9) 



with the Cj^o- given in [24]. Here, the renormalization is somewhat more tricky. 
It can most efficiently performed in a basis of a linearly independent subset of 
the operators dOf^ /dniq^ q = {u,d,s), as detailed in Ref.[24]. A good check 
on the rather lengthy expressions for the nucleon mass and cr7rAr(0) is given by 
the Feynman-Hellmann theorem, rh{drnN / dm) = O"7r7v(0), with m the average 
light quark mass, m = (m„ + md)/2. We remark here that in contrast to the 
order q^ calculation, the shift to the Cheng-Dashen point, (T7r7v(2M^) — (T7r7v(0), 
is no longer finite, i.e. there appear t-dependent divergences. We therefore do 
not consider this cr-term shift in what follows. We will also not discuss in detail 
the two kaon-nucleon cr-terms, (Tkn v), for similar reasons in this letter. 



4 Resonance saturation 



Clearly, we are not able to fix all the low-energy constants appearing in 
from data, even if we would resort to large Nc arguments. We will therefore use 
the principle of resonance saturation to estimate these constants. This works 
very accurately in the meson sector [25-27]. In the baryon case, one has to 
account for excitations of meson (i?) and baryon [N*) resonances. One writes 
down the effective Lagrangian with these resonances chirally coupled to the 
Goldstones and the baryon octet, calculates the pertinent Feynman diagrams 
to the process under consideration and, finally, lets the resonance masses go 
to infinity (with fixed ratios of coupling constants to masses). This generates 
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(a) (b) 

Fig. 1. Resonance saturation for the masses. In (a), a baryon resonance (from 
the decuplet) is excited, whereas (b) represents the t-channel meson (scalar or 
vector) excitation. Solid and dashed lines denote the octet baryons and Goldstones, 
respectively. 

higher order terms in the effective meson-baryon Lagrangian with coefficients 
expressed in terms of a few known resonance parameters. Symbolically, we can 
write 

C,s[U,B,R,N*]^jC,s[U,B] . (10) 



Here, there are two relevant contributions. One comes from the excitation of 
the spin-3/2 decuplet states and the second from t-channel scalar and vector 
meson excitations, cf. Fig.l. It is important to stress that for the resonance 
contribution to the baryon masses, one has to involve Goldstone boson loops. 
This is different from the normal situation like e.g. in form factors or scat- 
tering processes. Consider first the decuplet contribution. We treat these field 
relativistically and only at the last stage let the mass become very large. The 
pertinent interaction Lagrangian between the spin-3/2 fields (denoted by A), 
the baryon octet and the Goldstones reads (suppressing SU(3) indices) 

JCab^ = T e,^(Z) B + h.c. , (11) 

V 2-^p 



with T the standard ^ — > | isospin transition operator, Fp — 100 MeV the 
(average) pscudcoscalar decay constant and C = 1.5 determined from the 
decays A Bit. The Dirac matrix operator G^,/(Z) is given by 



(Z) = ^^,-(z+ 1)7^7, . (12) 



For the off-shell parameter Z, we use Z — —0.3 from the determination of 
the A contribution to the nN scattering volume 033 [28]. The tadpole graphs 
shown in Fig. lb with an intermediate vector meson only contribute at order 
q^. This is evident in the conventional vector field formulation. In the tensor 
field formulation, the vertex Tr{Vni,u^u'^) seems to lead to a contribution at 
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order q^. However, in a tadpole graph one needs the index contraction n = v 
and thus has no term at order q^. Matters are different for the scalars. Denoting 
by S and 5*1 the scalar octet and singlet with Ms — Ms^ — 1 GeV, respectively, 
the lowest order coupling to the baryon octet reads 

Cfs = Ds Ti{B{S, B}) + Fs Ti{B[S, B]) + Dg, S, Tr{BB) (13) 



where the coupling constants Ds, Fs and Ds^ are of the order one. ^ In fact, 
there are no empirical data to really pin them down. From the decay pat- 
tern of the low-lying baryons one can, however, estimate these numbers to be 
small. We will generously vary them between zero and one. For the couplings 
of the scalars to the Goldstones, we use the notation of [25] and the parame- 
ters determined therein. Putting pieces together, all low-energy constants are 
expressed via resonance parameters and the baryon masses take the form 



L^^, 5 = e = ^(2Z2 + Z-l) , (14) 

with rn^B the contribution of 0{'w?J'^) and — 1.455 GeV is the average 
decuplet mass. Notice that it is convenient to lump the 0{mq) and the C(m^) 
corrections togther as it was done in Eq.(14) (in the and -Df). We have 
kept explict the baryon mass in the chiral limit. Of course, in the fourth order 
terms it could be substituted by the corresponding physical values. At second 
order, however, we would get a state-dependent shift, see e.g. [29], and we thus 
prefer to work with m. Alternative representations for the baryon masses are 
given in [24]. Let us briefly explain the origin of the various terms in Eq.(14). 
The A contributions are tadpoles with 1/m insertions from C^^b- Similarly, 
the f3 and e terms stem from tadpole graphs with insertions proportional to 
the low-energy constants 6o,i:»,F and 6j, respectively. Note, however, that in the 
resonance exchange approximation not all of the ten bi are contributing. The S 
terms subsume the contributions from these are proportional to the low- 
energy constants di. Finally, the terms of the type are the scalar meson 
contributions to the mass. They amount to a constant, state-dependent shift. 
These consist of terms of the types ~ Mj>, ~ Mf> In M|, and ~ Mj>M^ lnM|, 
see [24]. To be specific, we give the coefficients D'^'^''' and A for the nucleon. 



-l{D'-2DF + 3F')Mt.ln^y 



In what follows, we neglect the singlet field. 
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M2 
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_ FfM'^K + — (223D^ - 318DF + 351F2)M^ 

4 Qfi 
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415 



288 



96 
83 
72' 



72 ^ 



(15) 



The corresponding coefficients for the A, S and S and also the can be 
found in [24]. The tree contribution from C^^^^^^^^^ is subsumed in the D^. 

The numerical values of the As, D^,D^^, D% and are given in table 1 
(using F = 0.5, D = 0.75 and from [25]). We see that the dominant terms 



B 


\b [GeV^] 




D% [GcV2] 


D% [GeV2] 


D| [GeV] 




0.0049 


-39.875 


-2.3617 


0.0321 


-0.017 L>5 + 0.061 F5 


A 


0.0179 


-68.044 


1.0831 


0.0617 


-0.047 L>5 - 0.010 F5 


S 


0.0144 


-132.90 


-11.273 


0.1422 


+0.051 + 0.010 




0.0216 


-126.46 


-5.4799 


0.1324 


-O.O37D5- 0.082^5 



Table 1 

Numerical values of the state-dependent coefficients in Eq.(14). The are dimen- 
sionless. The Do are for C(i, from [25] and M5 = 1 GeV. 



at 0{m^) are indeed the tadpole graphs with an insertion from C^^j^^^^^^ (this 
holds for the masses but not for (7^Ar(0)). It is also instructive to compare the 
values we find from resonance exchange with the ones previously determined 
from KN scattering data. We have transformed the results of Ref . [23] into our 
notation. As can be seen from table 2, most (but not all) coefficients agree in 
sign and magnitude. Note, however, that this is only a subset of the coefficients 
considered in this work. The full list will be given in [24] . We remark that the 
procedure used in [23] involves the summation of arbitrary high orders via a 
Lippmann-Schwinger equation and is thus affiicted with some uncertainty not 
controled within CHPT. 
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J, 

t>i 




I, 

03 




oo 


I, 

bp 


od 


Reso. {Ds = Fs 


= 0) 


0.084 


-0.144 


0.108 


-0.216 


-0.738 


-0.264 


0.317 


Reso. {Ds = Fs 


= 1) 


0.100 


-0.111 


0.125 


-0.239 


-0.733 


-0.208 


0.345 


Ref.[23] 




0.044 


-0.145 


-0.054 


-0.165 


-0.493 


-0.213 


0.066 



Table 2 

Some low-energy constants from in GeV~^. In the first row, only the decuplet 
contribution is given. In the second row, scalar meson exchange is added. 

5 Results and discussion 



The only free parameter in the formula for the baryon masses, Eq.(14), is the 
mass of the baryons in the chiral limit, m, since all low-energy constants are 
fixed in terms of resonance parameters. In particular, in contrast to Ref. [7], 
the parameters bo, bp, and bp are no longer free. Also, at quadratic order in the 
quark masses the ambiguity between m and 60 is resolved, it is not necessary 
to involve any one of the cr-terms in the fitting procedure [5,7]. In fact, one can 
not find one single value of m to fit all four octet masses, mjy = 0.9389 GeV, 
mA = 1.1156 GeV, = 1.1931 GeV and = 1.3181 GeV exactly. We 
therefore fit these masses individually and average the corresponding values 
for m. The contribution from scalar meson exchange only enters the uncer- 
tainties of the numbers given. This is justified since the numerical values for 
the couplings Fs and Ds are supposedly small. A more thorough discussion 
on this point can be found in [24]. We have m^= 711 MeV, m/^= 679 MeV, 
ms= 877 MeV, m^= 728 MeV, with the following average 

m= 749 ± 125 MeV . (16) 

This number is compatible with the one found in the analysis of the pion- 
nucleon cr-term, where approximately 130 MeV to the nucleon mass were at- 
tributed to the strange matrix element ms{p\ss\p) (with a sizeable uncertainty) 
[30]. The spread of the various values is a good measure of the uncertainties 
related to this complete calculation. Let us take a closer look at the quark 
mass expansion of the nucleon mass, in the notation of Eq.(l), 

= 711 + 202 - 272 + 298 MeV = 939 MeV . (17) 

This looks similar for the other octet baryons. We conclude that the quark 
mass corrections of order nig, w?J'^ and are all of the same size.Q In Ref. [16], 

^ Note that for certain values of the scalar couplings, the fourth order term ro'^ 
can be significantly smaller. However, the nucleon mass is much more sensitive to 
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it was argued that only the leading non-analytic corrections (LNAC) ~ m^/^ 
are large and that further terms like the ones ~ are modestly small, of the 
order of 100 MeV. This would amount to an expansion in M|-/(47rFp)^ ~ 1/6 
with a large leading term. This expectation is not borne out by our results, 
the next corrections are as large as the LNACs. These findings agree with 
the meson cloud model calculation of Gasser [3]. A last remark about the 
baryon masses concerns the deviation from the Gell-Mann-Okubo relation, 
^GMO = (3mA + ms — 2miv — 2ms)/4 which empirically is about 6.5 MeV. We 
find Aqmo = 31 MeV, which is larger in magnitude than the value found in 
[5,6]. We remark that in our case the decuplet contribution is contained in the 
contributions and not in the m^/^ as in [5] . Therefore, in our case, Aqi^q is 
dominated by the piece. The sizeable uncertainty in the chiral limit masses, 
Eq.(16), does not allow for a very accurate statement about this very small 
quantity. It is also very sensitive to the scalar couplings. To get a better handle 
on this issue, one either has to be able to fix all pertinent low-energy constants 
at order from data or improve upon the resonance saturation estimate used 
here by including e.g. the mass splitting within the decuplet and the SU(3) 
breaking of the decuplet-octct-mcson couplings. A better understanding of this 
topic is, of course, at the heart of the determination of the quark mass ratio 
R — {rrid — rriy) / (m^ — m) from the baryon masses (once the electromagnetic 
corrections have been included). 

The pion-nucleon ci-term is completely fixed. Using for m its average, we find 
(no scalar resonance contribution) 

(7^^(0) = 48 MeV , (18) 

with an uncertainty of about ±10 MeV due to the spread in m. This number 

compares favourably with the one found in Ref. [30], (Tt^n{0) = 44 ± 9 MeV. 
We stress that this result reflects a very non-trivial consistency for the com- 
plete calculation to quadratic order in the quark masses using the resonance 
saturation principle in the scalar sector. The additional contribution from the 
scalar meson exchange is accounted for in the ±10 MeV uncertainty. To be 
specific, we have cr^Af(O) = (48 + Ds ■ 0.5 — Fs ■ 2.0) MeV. It is furthermore 
instructive to disentangle the various contributions to (7^Ar(0) of order q^, 
and g^, respectively, 

(7^^ (0) = 54- 33 + 27 MeV = 48 MeV , (19) 

which shows a moderate convergence, i.e. the terms of increasing order become 
successively smaller. Still, the contribution is important. Also, at this order 
no TTTT rescattering effects are included. We notice that using the values for 

the scalar contribution than the other octet masses. 
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bo.D,F and bi as determined in Ref.[23] leads to a much increased fourth order 
contribution. 



6 Summary and outlook 

In this paper, we have used heavy baryon chiral perturbation theory to cal- 
culate the octet baryon masses to quadratic order in the quark masses, in- 
cluding 20 local operators with unknown coeffcients. These low-energy con- 
stants were fixed by resonance exchange. The dominant contribution comes 
indeed from the excitation of the spin-3/2 decuplct fields. Tadpole graphs with 
scalar meson exchange only lead to small corrections. This left us with one 
free parameter, the baryon mass in the chiral limit, which could be determined 
within 18% accuracy, and is compatible with the value inferred for the strange 
matrix-element msiplsslp) in Ref.[30]. Furthermore, the pion-nucleon cr-term 
comes out surprisingly close to its empirical value, cr^Ar(O) = 48 MeV, with an 
uncertainty of about ±10 MeV. This first exploratory study of the three 
fiavor scalar sector of baryon CHPT points towards a significant improvement 
compared to previous investigations which were mostly confined to so-called 
"computable" corrections and/or fitted a few of the pertinent low-energy con- 
stants. However, the calculation is not yet accurate enough to determine the 
quark mass ratio R reliably from the octet masses. Furthermore, we did not 
address the kaon-micleon a-tcrms, the corresponding shifts to the pertinent 
Cheng-Dashcn points together with the strangeness content of the proton here. 
We will come back to these topics in Ref . [24] . 
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